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ABSTRACT 


The effects of wind-shears on the approach performance of 
a STOL aircraft are ana Iv red usinq the optimal-control model of 
the human operator. This analysis involves a time-varyinc 
situation that is more complex than is traditionally treated 
by human operator modellinq techninuef. The extensions to 
the time-varvinq case are discussed aid results are presented 
that Illustrate the effects of wind shears on cateqory II 
window-performance and the effects 0‘i performance of pilot 
tinc-delav and of variations in pi let qain durinq approach. 

1. IJITRODUCTION 

In an instrument approach-to*landing the pilot attempts 
to maintain his position on the qJide-path in the presence of 
external disturbances. An important part of the task involves 
compensation for errors introduevd by winds. Most analytic 
studies of the problem based on human operator models have 
considered only the effects of sero-mean, random turbulence. 
However, winds with a non-zero mean component can be quite 
significant and may provide the dominant problem. These mean- 
winds" generally vary with altitude. The rate of change of wind 
speed with altitude is called the shear variation and the alti- 
tude dependent winds are often refered to as wind-shears. The 
mean-wind can be described in statistical terms, i,e,, the wind 
direction and speed are, in general, random variables. However, 
in a given approach-to-landinq, a specific ’•sasple" mean-wind 
is encountered. It is the response to particular samples that 
is usually of interest, rather than the response to the distri- 
bution as a whole. 


'This research was performed under Contract NA32-<;652 for the 
NASA, Ames Research Center, 
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The introduction of wind shearn into the approach problem 
complicates the sltuat><on considerably from the standpoint of 
analysis by human operator modelling techniques. The distur- 
bance input no longs:: has stationary statistics and, furthermore, 
one expects the pilot to be capable of some "prediction" or 
other high level adaptation with respect to the shear variation. 
These two factors tend to compromise the theoretical and experi- 
mental bases for most quasilinear describing function models 
(1), although interesting attempts have been made to extend 
those models to ouasi-prediotable inputs (2), 

Here, the optimal control model of the human operator (3) 
is used to analyse approach in a wind- shear environment. It is 
relatively straightforward, theoretically, to apply the optimal 
control model to this problem, because of the model*s time domain 
foundation and its normative nature. This was demonstrated in a 
previous application to time-varyinc approaches involving constant 
updrafts [4], That study also showed that model results were 
reliable predictors of experimental data for the time- varying 
situations analysed. 

This paper is an extension rf the analysis presented in (4); 
specifically, a more general clast of disturbances is considered, 
Mi some time-varyinc facets of the optimal control model are 
examined in greater Mtail, The results are believed to be of 
interest both in terms of predicting the effects of wind-shears 
on approach performance and in leading to further understanding 
of human operator modelling techniques based on optimal control. 

On the other hand, unlike experimental data to confirm or 
reject the model results are, fortunately, not available at 
this time. 

The paper begins with a brief discussion of the modifica- 
tions to the human operator model required for analysis of the 
response to wind shears. Then, results are presented for a 
STOL approach that illustratei U> the effects of wind shears 
on performance at the approach "window" ^ (2) the effect of 
pilot time-delay performance ; and (3) the effects of changing 
pilot gain during approach, A more detailed presentation of the 
ideas and results may be found in (S), 

2, MODIFICATIONS OF OPTIMAL CONTROL MODEL FOR TlMB-VARYIir, 

DISTURBANCES 

Consider the lin*»dr dynamic system 

x*Ax*Bu4-Ew-Fs (1) 


•419-. 
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where x is a vector representing the system state, n is a vector 
of control inputs, w is a vector of rero-mean, gauss'ian, white 
noises , 2 a vector o7 time-varying input disturbances and A, B, 

E and r constant matrices of appropriate dimension* We assume 
that z satisfies 

z « A 2 ; z (t^) “ 2 ^ (2) 

with A a constant matrix, (Note that if B » F then z can be 
used to model a mean component of w*) Varrous~disturFanee 
inputs may be oenerated from the model of Equation (2) , parti* 
cularlv if impulses or lumps in **initial ^'^nditions” are 
allowed ( 6] • 

To examine how the control structure of the optimal control 
model 18 modified to account for the disturbance (2)# we iqnore 
the limitations on the human's observation processas, (time-delay 
and observation noise) and assume the As and z(to) are known. 

This implies the disturbance-state, 2 (F) , is Fnown for all t. 

This assumption is unrealistic and is made here only for exposi- 
tory reasons; in applying the model, £\t) is estimated from the 
available noisv, delayed data, as are'^the other system states. 

The operator's control u is assumed to be chosen to mini- 
mize the ouadratic cost functTonal 

J i'2 * ♦ 


Where 0 and R are positive, semi-definite matrices and G is 
positive defTnite, The solution to this problem is welT-known 
and IS developed in detail in the control literature [6-8] • In 
particular, it can be shown that the optimal control law satis- 
f j es 


Tj^(t)u* ♦ u* - - ?l<t)x K(t)z ) (4) 


Where and K are time-varyinq matrices that are obtained 

bv sol^no an appropriate matrix Ricatti equation (see (8) and 
[3] for details). Moreover, if the input disturbance does not 
have an exponentially qrowinq component (i,e,, if the eigen- 
values of A are < 0) , the control gains in (4) approach a 
constant as T>,«. ” Thus, when the system matrices and cost 
functional weightings are time-invariant and when T is large 
relative to system time-constants, we may take advantage of 
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the enormous eimpli fleet ion afforded by a constant gain solution,* 
Under these conditions Equation (4) mey be written as 


U* • - (Tj^ 8 ♦ X) <£• JC 4 K 2) (S) 


whers ^ and L* are identical to the neuro-motor lag and optimal 
faedbaek gain*matrieaa obtained for the optimal control model of 
the human operator under atationary oonditiona (31, Thus, the 
modifioation to the baaie medal to account for the tima-varyinq, 
mean dieturbenoe is a aet of feedforward gains aotino on the 
(estimate of) the diaturbanca,** This raault, which is a direct 
eonaaquanea of the normative aaeumption, is satisfying intuiti- 
vely and is not ineenaiatant with the "Sueeeaaive Organization 
of Paroaption" concepts diaeusaad in, a,g,, [9), 

As mentioned previously, z(t) is to bo aatimatad from the 
available ebaervatien data, Xn*tarma of the estimation problem, 
it will be aaaumed hare that the initial condition £0 ia a oauesia 
random variable with known mean and variance. Thus, sit) is 
gauasian and its distribution is known, for ail t, T^a actual 
value of 2 (t) cerraaponding to a specific sample disturbance from 
the distribution in not known and is aatimatad from the displayed 
variablaa by maana of a Kalman filter, 2t ean be shown (S) that 
if the aampla path does not eorraapond to the path generated by 
the mean of Sp, then the error and the aatimata are correlated, 
the filter being optimal only in terms of the distribution of 
Sq, The expressions used to compute the corresponding estimation 
errors for this ease are <*evalored in (5), 

Finally, it is important to note that the man-machine sys- 
tem response to a sample disturbanca z(») is a random variable, 
even when w 2 0. The reason for >.hiB**is the human operator's 
randomnasa”whioh is raflaetod in the observation and motor 
noises of the optimal control model. 


•Actually, the ^solution* may not be a eolution at all, in that 
the infinite-time coat may be unbounded, Kewavar, this is a 
minor ineonvanianea that ean be eireumvanted readily xn 
applying the reaulta (S, 7), 

••An altarnativa, perhaps lass aatiafvine, interpretation is 
that the control law involvaa integral as wall as proportional 
faedbaek (6, 8| • 
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ANALYSIS OF APPROACH PERFORMANCE 
3,1 Approach Scenario 

We consider the longitudinal approach of an Augmentor Wing 
Jet STOL Research Aircraft (AWJSRA) , the C*8M.* The aircraft is 
assumed to be initially on the nominal n«S deq.) qlide*8lope 
With a nominal airspeed of approximately 31 m/s. Linearized 
perturbation eauations were used to describe the aircraft motion. 
The pilot was assumed to control elevator and * nozzle" (thrust 
vector) in a continuous manner^ whereas throttle was assumed to 
be fixed at the appropriate trim setting. The corresponding 
equations of motion^ in state variable formi are given in (SI, 

The basic scenario for the analysis of longitudinal approach 
performance involved starting at an initial range of 1500m with a 
constant wind velocity corresponding to the value h ® 125m 
(R z 1160r), In all cases^ turbulence having the Dryden spec- 
tral form (10) was assumed present. The scale-lengths of the 
turbulence were not varied with altitude; they were set at the 
constant value apnropriate to the decision height. Gust inten- 
sities corresponded to a value that would not be encountered 
more than 10% of the time, i.e,, a fairly severe condition. 

The pilot's disnlav was the EADl status display used in 
the STOLAND proorair (ill and described In another paper presented 
at this conference (12) • The STOLAND-status display shifts from 
an annular presentation of glide-path error to a height presen- 
tation at Range z S75n (h z 75m z 250 ft,). The model for the 
human operator takes this variation in display gain into account 
bv modifying appropriate observation parameters. Values for the 
parameters of human operator model were the same as those used 
for a corresponding steady-state analysis of performance at the 
decision height and are given in (5) and (12), Suffice to say, 
that the values for parameters corresponding to human limitations 
were essentially the same as in previous studies (3, 4) and task- 
related parameters were determined from analysis of this problem 
( 5 , 12 ). 

In addition to the above "basic" condition, the effects of 
variation in pilot gain were analyzed. It had been noted in (4) 
that the pilot might vary his gains so as to tighten control as 
the decision-height was approached. An analogous result can be 
obtained with the model by making the cost functional weightings of 
Equation (3) range-dependent. Indeed, if it is assumed tliat the 


•An analysis of the lateral performance is given in (5), 


pilot attaches a fixed penalty to angular deviatione from the 
glide-path, rather than to linear deviations, then the weioh- 
tlng on height-errors will be range-dependent. We investiisted 
two conditions with respect to this "gain-sehedulino"; 

1* Constant Gains in which the cost functional 
weightings were constant and corresponded 
to Category XI "window" performance require- 
ments at the decision height, 

2, Varying Gains in which angular glide-path 
error toierancas were assumed to be constant 
(corresponding to the allowable Category XI 
window error) , For this case, weightings 
(and, hence, gains) were changed in three 
stages according to the following schedule; 

1500 < R < 1160m i • ,0029, * *0326 

1160m ? R ? 575m t gj? ■ *0117, gP • ,133 

575m < R 2 230m i gJJ • ,073, qg ■ .83 

Thus, the weighting over a range-interval corresponded to the 
weighting appropriate to the end-point of that interval, a con- 
servative choice* The intervals were chosen, as a matter of 
convenience, so that the end-points corresponded to points 
where other changes in the approach scenario were reouired* 

A limited examination of the offset of pilot time delay 
was also conducted* It was expected that t)*e human's time 
delay would increase scores but would not liter the basic 
character of the results* Because inclusion of the time delay 
increases significantly the costs of the time-varying computa- 
tion, we decided to assume the time delay was zero* However, 
a comparison case was obtained with a time delay of *2 see, 

(the nominal value found in previous studies) to illustrate 
the differences one might expect from including time delay. 


3*2 Modelling the Wind-Shears 

‘ In modelling wind-shears both dynamic and kinematic effects 
should be considered. Xn addition, from the standpoint of the 
optimal control model, ii iB desirable (though not necessary) to 
convert the altitude dependence of the shear to an equivalent 
tire dependence. In this section, kinematic effects of wind 
shears are modelled as is the conversion to a time-dependent 
wind; dynamic effects are accounted for by considering the 
shear-components in the same fashion as turbulence velocity 
componente (5), Only horizontal wind-shears are considered. 
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Figure 1 illustrates the pertinent geometry# The aireraft*e 
altitude (h) is given by 


h » ♦ <5h » K tan + (Sh (6) 


where h„ is the "nofivinal* altitude, i#e., the altitude of the 
glide-slope at the aircraft's range, H, and 6h is the altitude 
error# The rate of change of the nominal altitude may be 
expressed in terms of the ground speed <or range-rate) • 

h. ■ R tan r a (U^ u) tan F (7) 

n o o o 


where 

u » x-bodv axis component of perturbation in ground speed 

» x-body axis component of nominal airspeed 

The aircraft's sink-rate is 

h » V sin y » ♦ u) sin 4- hy) 

o o 

« (t’ + u) (sin F ♦ COB r • Ay) 

o o o 

and 

* • • 

6h - h - h^ 
n 

* (IJ + u) cos r • Ay • u) cos r • (6-e) 

o o o o 

or 

5h ® (U^ ♦ u) cos • 0 - - 9^ — H cos • w 


Cnvxtion (8) is used to account for kinematic effects of shears# 
Kov/ever, this enuation is nonlinear and the product terms ^ U0 
anc uw, mav not be negligible if u is a significant fraction of V • 
7!ixp IS the case for ti«e v^iiu's to be considered here# To main- ® 
tain linearity and reduce the errors associated with neglecting 
the product terms, the "average" wind-velocitv during approach 
was substituted for u in Eouation (8) i this tende to minimise the 
maximum error associated with assuming constant ground speed# 

The wind-shears to be considered here are enumerated in 
T^ble 1# These winds are idealisations of more exact models for 
mean-winds# They were us i in this analysis so as to be compatible 
with a concurrent simulation study# he now show how these winds 
mav be represented as time-varying disturbances r it turns out that 
this can be done with considerable fidelity# 
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+ Indicates tailwind or crosawind frcan left side* 
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• 

This equation alone with the dvnamical equation for u allows the 
wind shear to be expressed in terms of other^ non*input related 
state variables* An even simpler representation is possible* 

If the pilot is maintaininq airspeed reasonably well^ then 


u 


^MW ®o 


( 12 ) 


and 9 the state-variable representation for the mean wind distur- 
bance is given by 


«1 “ «MW “ *2 > *1^®> “ 


B (a tan cos 0^)22 ^ *2 “ * ^o ^^o 

®o 

where Z 2 ^ is the mean-wind velocity at the onset of the shear* 

3*3 Results 

Mean and standard deviation scoresi at the decision-heieht 
(31m) ^ are compared in Table 2* Several points are worth mentionino 
First# the constant-gain no shear# zero-delav results are virtually 
equal to those of a corresponding steady-state analysis (5) • This 
is more than a chmek on the computer programi it shows that in the 
absence of shears# the approximately 12 50-1 300m approach distance 
is sufficient for the errors to reach steady-state* Second* the 
effect of the wind-shear is more than just a non-zero mean res- 
ponse* It may be seen that the standard deviation of the tracking 
errors and of the controls is increased* This is a result of the 
coupling in the model of mean-and variance-responses that arises 
from the dependence of the observation noises on the rms signal 
values and that of the motor-noise on rms control* in terms of 
missed approach probabilities# the increase in variance is the 
more significant effect* Third# the effect of time delay is# as 
expected# to increase mean and standard deviation of the error* 

The magnitude of the effect is largest for heicht-error with 
approximately a 35% increase in mean and a 10% increase in 
standard deviation* 



-744 




Table 2 


PERFORMANCE AT DECISION-HEIGHT FOR 
VARIOUS ANALYSIS CONDITICmS 



The final effect illustrated in Table 2 is that resulting 
from allowina the gains to varf. When compared with the constant 
gain case, it is seen that the principal effect at the window 
js on the mean-response . Thir effect, though large percentage- 
vise, i;: v’rtually necrliaiblc in terms of the missed-approach 
rrobabilitv. The differences between constant and varying gains 
are more p* >nounced in the "tiRie-histories** shown in Figure* 2t 
These time-histories are curves passed through data points 
obtained every 5Cm in range. The jump-discontinuities for the 
varyinq-qain case arise from the instantaneous gain-change and 
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FIGURE 2, Effect of Varyi .1 Gains on Approach Trajectories 
(a) Height and Nozsle Response 
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FIGURE 2. Effect of Varying Gains on Approach Trajectories 
(b) Pitchr Airspeed and Elevator Response 
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Ae associated jump in control value. These "jumps" a^arently 
iecay very rapidly* Because heicht errors are weighted less^ 
at i»re distant ranges, they are allowed to build up to a greater 
extent in the varying-gain case; however, as the threshold is 
approached the errors begin to be reduced rapidly (because of the 
higher weighting) , so that window perfonnance is not significantly 
different for the varying- and constant-oain cases (Table 2). 

Apart from differences in height control, the principal difference 
between the two-cases is in the initial transient in elevator and 
pitch. It se«ns clear that the early reduction in height errors 
for the constant-gain case is a result of a rapid pitch-down. 

The excellent "window" perfonnance obtained in the above 
analyses is somewhat misleading* As can be seen from Figure 2, 
the nozsle limit (29*) is less than one standard deviation from 
the mean for much of the approach (after the wind-velocity starts 
changing). Thus, a high percentage of the time the nozzle will 
exceed its limit* What this means is that the rate of <3es^nt 
capability of the aircraft, with throttle fixed , is insufficient 
for this wind* Further, the wind is of sufficient severity ^ 
place the entire linearized analysis in ouestlon* On the other 
hand, the analysis suggests that suitably scaled-down winds 
may be adeouately controlled by nozzle and elevator inputs alone* 

In an attempt to get some estimate of the control-limited 
performance for the decreasing tailwind, a trajectory was obtained 
for a case in %diich nozzle control and control-rate were heavily 
penalized in the region where excessive nozzle-control had been 
observed, i.e*, in the ll60m < R < 575m Interval. (Weightings 
on nozzle and nozzle-rate were multiplied by 50), To 
transient effects resulting from the initial constant wind to 
die out, the approach was started at 2000m* The result for 
height-error and nozzle-position is shown in Figure 3. It can 
be seen that nozzle responses to the shear variation in the 
heavily penalized region are virtually nil and the height errors 
increase aecordinaly. When the penalty is reduced, R < 575i», a 
relatively large mean-nozzle motion ensues in an attempt to reduce 
the mean-error* While some reduction occurs, the mean height 
error at the decision height is still three times the allowable 
Category II error. Although these results arc not intended to 
be definitive, they do illustrate the problem posed by this wind, 
when the throttle is fixed. 

A constant-gain trajectory for the increasing tailwind was 
also obtained and the results are shown in Figure 4.* windw 

performance for this wind is compared with that for the decreasing 
tailwind in Table 3. Note that the turbulence intensity and 


Xs can b« Men (di>e"^nd appro*im.tion la not aa cloaa to th. 
Idoalited wind aa for the previous casf , but certainly qood 
enough. 
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FIGURE 3 
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FIGURE 4. Response History for Increasinq Tailwind Shear 
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spectrum IS t.he same for the two cases. It mav be seen that 
heiQht errors are controlled more effectively for the increasino 
tailwind; airspeed is less-well controlled. The overall effect 
is a definite improvement, as could be expected. Two other 
points are worth notinc* Referrinc to Table 3, we find that the 
standard deviation of the heiaht and sink-rate errors for the 
increasino tailwind are very close to those obtained for the 
steadv-state (no-delav) case. Thus, it appears that with the 
increasi.'.g tailwind (which starts out as a headwind), enough 
tiifte 19 added to the aoproach to allow the "pilot" (model) to 
reduce the errors to values comnensurate with an approach of 
infinite lenoth. The second point is that the noasle-control 
reuuirements are not so excessive (in relation to capability) 
as for the decreasing tailwind. Thus, one might expect these 
1 ^sults to correspond more closely to a realistic situation. 


S . CONCLUS rON 

Some aspects of the STOL approach in a mean-wind with 
snear-variation c^^^e been analvzed with the optimal control model 
of the hrman operator. Results were obtained for longitudinal 
control with an EADZ status display. In general, the wind-shears 
dearade performance by producing both mean errors and increased 
variability in the response, with the increased variability 
appearing to be the major effect. 

Two wind profiles were considered in the analysis of longi- 
tudinal control, a decreasing and an increasino tailwind. Rela- 
tively good performance at the window was obtained in both cases, 
however, for the decreasing tailwind, the results showed that 
with the throttle fixed, excessive nozzle-control was reauired 
for wind compensation. When the nozzle control was limited 
(indirectly, bv penalizing control motions subseauent to shear- 
onset), the height errors increased significantly. The relati- 
vely good performance for the increasing tailwind was achieved 
with control reouirerents that were not so excessive and, conse- 
auentlv, represent a more reliable result. The better performance 
IS undoubted!^ due to the additional time available for error 
compensation and is, of course, to be expected. 

To the extent that thev were investigated, the wind-shear 
responses tended to confirm essentially the results of a corres- 
ponding steady-state analysis (5] (albeit that performance was 
worse in the wind-shear). The details of the transient responses, 
may be easily explained. Thev depend very much on the specific 
assumptions about initial conditions and on pilot strategy, which 
IS not the least hit surprising. If one is interested in 
cino or predicting a particular time history (ensemble) then it is 
essential that conditions used in the model match those of the 
experiment. 
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With respect to 'uture %raz)c, it appears most important to 
validate the results of this analysis with simulation data* 
Approach data is needed to pin-down details of the pilot's time- 
varying adaptation* Transient data for a sufficient number of 
runs to provide reliable statistical information would be most 
helpful* 


Table 3 

COMPARISON OF WINDOW PERFORMANCE FOR 
DIFFERENT TAILWINDS 
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